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Abstract 

In this work, the controversy between the interpretations of recent measurements on dense 
aluminum plasma created with Linac coherent light sources (LCLS) X-ray free electron laser 
(FEL) and Orion laser has been addressed. In both kind of experiments, helium-like and 
hydrogen-like spectral lines are used for plasma diagnostics . However, there exist no precise 
theoretical calculations for He-like ions within dense plasma environment. The strong need 
for an accurate theoretical estimates for spectral properties of He-like ions in strongly coupled 
plasma environment leads us to perform ab initio calculations in the framework of Rayleigh-Ritz 
variation principle in Hylleraas coordinates where ion-sphere potential is used. An approach to 
resolve the long-drawn problem of numerical instability for evaluating two-electron integrals with 
extended basis inside a hnite domain is presented here. The present values of electron densities 
corresponding to disappearance of different spectral lines obtained within the framework of 
ion-sphere potential show excellent agreement with Orion laser experiments in Al plasma and 
with recent theories. Moreover, this method is extended to predict the critical plasma densities 
at which the spectral lines of H-like and He-like carbon and argon ions disappear. Incidental 
degeneracy and level-crossing phenomena are being reported for the hrst time for two-electron 
ions embedded in strongly coupled plasma. Thermodynamic pressure experienced by the ions 
in their respective ground states inside the ion-spheres are also reported. 


1 Introduction 

The study of confined quantum mechanical systems has attracted immense attention from researchers 
around the world due to the novel and unusual structural properties exhibited by such systems when 
subject to spatial limitation [1]. A wide variety of physical situations are manifested in nature 
that relates to spatially confined systems such as atoms or molecules trapped in zeolite sieves [2], 
fullerenes [3], plasma environment |1], solvent environment [S], under high pressure in the walls of 
nuclear reactors [6], quantum dot or artificial atom [^, molecular containers, storage of fuel cells 
matter under high pressure in Zovian planets [lO] etc. Along with the experimental and technolog¬ 
ical development, theoretical research plays a fundamental role for designating appropriate models 
in order to explore and predict the behavioral changes of a conhned system. The present study 
is focused on atomic systems embedded in plasma environment. In recent years, atoms placed in 
external plasma environment have received considerable attention from researchers [TT1 - [T9] due to 
their wide applications in various disciplines of science e.g. astrophysics, condensed matter physics, 
biology etc. While dealing with plasma that follows classical statistics, a coupling parameter (T) 
defined as the ratio of the average electrostatic energy and the average thermal energy is introduced. 
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r < 1 corresponds to weakly coupled plasma (WCP) for which the effective potential experienced by 
the embedded ion is expressed according to Debye model [SU] and F > 1 denotes strongly coupled 
plasma (SCP) where the potential is taken from Ion-sphere (IS) model [21]. According to the IS 
model, a sphere (termed as Wigner-Seitz sphere) surrounding a positively charged ion is considered 
in such a way that the plasma electrons within the sphere neutralize the positive ion. The size of 
the Wigner-Seitz sphere will decrease when the number density of plasma electrons (rig) increases. 
The temperature (T) of the plasma does not appear directly in this model but it is implicit as n^. 
is different for different temperatures. The domain of the effective potential representing the SCP 
surrounding is hnite in case of IS model in contrast to the long range character of the screened 
Coulomb potential used in Debye model ra for WCP environment. The examples of WCP’s are the 
gaseous discharge plasma {T ~ 10^ K and n ~ 10^^/c.c), plasma in controlled thermo-nuclear reac¬ 
tion (T ~ 10® K and n ~ 10^®/c.c), solar coronal plasma (T ~ 10® — 10® K and n ~ 10® — 10^®/c.c), 
Tokamak plasma {T ~ 10® — 10^ K and n ~ 10® — 10^®/c.c) etc. SCP’s (temperature varies and 
typical densities > 10^®/c.c.) are observed in highly evolved stars in high density states, interior of 
Jovian planets, explosive shock tubes, two-dimensional states of electrons trapped in surface states 
of liquid helium, laser produced plasmas etc. Spectral line shifts, pressure ionization, ionization 
potential depression (IPD) and line merging phenomena occur in SCP environment due to the defor¬ 
mation of the ionic potential by the plasma helds. Such properties and knowledge about ion-plasma 
interaction can effectively be utilized for diagnostics and the investigation of X-ray opacity of matter 
under conditions prevailing in stellar interiors. The experimental observations using laser produced 
plasmas for (7, Al and Ar by Nantel et. al. iza, Saemann et. al. [23] and Woolsey et. al. have 
explicitly demonstrated the effect of SCP on the spectral properties of such systems. The laboratory 
plasma conditions (T and rig) undergo rapid changes w.r.t. where local thermodynamic equilibrium 
is not maintained. Consequently, the experimental measurements become extremely complicated 
leading to a loss of accuracy and, till the end of the last century, this accuracy level was not even 
mentioned in most of the experiments. 

In recent years, a remarkable improvement has been made [2SH2B] with the advent of Linac Coherent 
Light Sources (LCLS) towards creation of relatively long-lived high-density plasma at homogeneous 
temperature and densities. In these experiments. X-ray free-electron Laser (PEL) was used to create 
plasma with densities up to almost one order higher than solid Al and then spectral line prohles 
of different charge states of Al were used for diagnostics. The effect of IPD on the emitted spectra 
as a function of is explored experimentally by observing the disappearance of spectral lines of 
H-like and He-like Al. During the observation of K-shell fluorescence of highly charged Al, Ciricosta 
et al. iza found that the IPDs measured were not consistent with the predictions of the most widely 
used theoretical model of Stewart and Pyatt (SP) [29] but in good agreement with an earlier model 
due to Ecker and Kroll (EK) [30]. However, this observation was questioned in a subsequent theo¬ 
retical study by Preston et al. [21] where detailed simulations were carried out for the spectral lines 
of H-like and He-like Al to study IPD by using both SP and EK (in a modihed form) models. In 
experiments, the intensities and Stark-broadened widths of He-/9 and Ly-/3 spectral lines are used 
for main diagnostics. A direct measurement of ionization potential depression is a difficult task be¬ 
cause of its indistinguishability from the effect of spectral line merging due to Stark broadening [32] . 
Hoarty et al. [331132] have been able to overcome this difficulty and their measurements for Al plasma 
using Orion laser are in closer agreement with SP model of IPD than the EK model. This situation 
clearly warrants extensive and accurate ab initio study of atomic structures within dense plasma 
environment. Very recently. Son et al. [35] have adopted a two step Hartree-Fock-Slater approach 
to assess the IPD effect for A/®’*' to Af~^ within plasma where a muffin-tin flat potential was used. 
The IPDs calculated by Son et al. [53] he between the SP and modihed EK models and in some 
cases, are close to SP model. But so far, no extensive theoretical calculation on IPDs for He-like 
ions has been performed. It should be noted here that both SP model and EK model for estimating 
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IPDs are derived within the framework of IS potential. The only theoretical work for He-like ions 
in the held of SCP by using IS potential is due to Sil et. al. [55] where both non-relativistic and 
relativistic calculations were carried out using time-dependent perturbation theory. Though Sil et. 
al. [36] demonstrated that the relativistic IS model yields consistent results in predicting the spectral 
line positions for the systems considered, some anomalies such as better agreement of non-relativistic 
results with experiments than relativistic ones are observed in their data [SB] . Such strange features 
may arise due to improper inclusion of electron correlation in basis set within a hnite region. A ma¬ 
jor challenge for precise theoretical calculations is, therefore, to develop an appropriate methodology 
where the effect of electron correlations within a hnite domain is aptly included. 

To the best of our knowledge, there exists no calculation of He-like atoms embedded in SCP us¬ 
ing Hylleraas type basis set though it is well accepted that within the framework of Ritz variational 
technique, explicitly correlated wave functions expanded in terms of Hylleraas basis (and its vari¬ 
ants) can produce most accurate non-relativistic energies of He-like atoms. These methods have been 
applied extensively to free He-like systems whereas for spatially conhned two-electron systems, such 
studies are limited to S states only [371H0] . According to Laughlin and Chu [39], the generalized 
Hylleraas basis sets used in such calculations suher the loss of linear independence for large dimen¬ 
sions of the wave functions and hence all the calculations [HTfHO] were limited to small dimensions 
(at best 25). Laughlin and Chu [39] made an ehort to remove this difficulty and extended the basis 
size up to 95 parameters where they have to compromise with the flexibility of the non-linear param¬ 
eters. Recently, for states of He-like systems under spherical conhnement, the present authors 
have calculated the energy values [H] by using standard Hylleraas basis set of dimension 161 and the 
results have been conhrmed by Montgomery and Pupyshev [12] • In the present work, a successful 
effort has been made to develop a general methodology in Hylleraas basis for both S and P states 
of He-like systems. The hnite domain two-electron integrals with hexible parameters are evaluated 
where the problem of linear dependency in larger dimensions is clearly avoided. 

We have estimated precise non-relativistic energy values of Isns (^S*®) [n = 1 — 3] and Isn'p {^P°) 
[n' = 2 — 4] states of He-like C, Al and Ar within SCP environment. Accuracy of the computed 
energy eigenvalues have been tested systematically over an extended range of parameters and also by 
increasing the number of terms {N) in the expanded basis sets. The plasma densities (ue) are varied 
from a low value that corresponds to almost a free system to a very high one that leads the ion to¬ 
wards destabilization. The plasma electron densities in different experimental conditions [22H2SIE3] 
are well covered within the density ranges studied here. The energy eigenvalues of ns (^S) [n = 1 — 2] 
and n'p (^P) [n' = 2 — 3] states of H-like C, Al and Ar in SCP are also estimated to determine the 
variation of ionization potential (IP) w.r.t. Ug- As Ug increases, both the two-electron excited states 
as well as the respective one-electron threshold move towards destabilization, thereby reducing the 
IP. It is remarkable that after a certain value of Ug, the two-electron energy levels move above the 
respective one-electron energy level and become quasi-bound. Incidental degeneracy [43] and subse¬ 
quent lev el-eras sing phenomenon between the excited states such as ls2s (^R®) and ls2p (^R°) under 
SCP have been observed. Such features are novel in the context of foreign atoms in SCP and being 
reported for the hrst time in literature. Due to spatial restriction imposed upon the wave function 
according to IS model under SCP environment, the ion will feel a pressure inside the Wigner-Seitz 
sphere. The variation of thermodynamic pressure w.r.t. plasma density is also calculated. The paper 
is organized as: an outline of the basic theory used and details on the evaluation of the basis integrals 
are given in Section 2, followed by a discussion on the results in Section 3 and hnally concluded in 
section 4 with a view towards further application of the present methodology in related helds. 
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2 Method 


The non-relativistic Hamiltonian (in a.u.) of a two-electron ion placed inside SCP environment can 
be written as 






(1) 


Visi^i) is the one-electron term of the modified potential energy as ^seen' by the i-th electron within 
plasma environment. It is to be noted that in this model, the electronic repulsion part in the potential 
is completely unaltered. The spherically symmetric potential Vis{ri) experienced by a positive charge 
ion surrounded by a one-component plasma within the ion-sphere [2T] is given by 


Visin) 


Z jZ-N,) 
n 2R 



( 2 ) 


where R is the Wigner-Seitz radius [21], i.e. the radius of the surrounding ion-sphere, Z is the 
nuclear charge and Nf> (< Z) is the number of electrons present in the ion. For helium-like ions, 
A^e = 2 is being taken. The Schrodinger equation i/T = is to be solved to obtain the energy 
eigenvalues where the wave function is subject to the normalization condition = 1 within the 

sphere. The structure of the potential demands that there is no electron current taking place through 
the boundary surface of Wigner-Seitz sphere, and the orbital wave function T satisfies the boundary 
condition 


T(r) = 0 at r > R 


(3) 


This boundary condition plays a significant role in behavioral changes of the confined atoms. The 
plasma electrons within the ion-sphere neutralize the central positive charge and the size of the 
Wigner-Seitz sphere is determined by the condition of overall charge neutrality that yields 

1 
3 


3(Z-iVe) 
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The above expression for ‘i?’ is used to determine the IPD according to the SP model [22]. However, 
in EK model [30] for determining the IPD, this radius was calculated in a somewhat different way 
where both the electron density (ng) and ion density (n,) are considered. According to EK model [30] . 
the radius of the sphere would be expressed as 


Rek — 
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(5) 


The EK model [00] is relevant when the ion density is appreciably high and can affect the mean 
separation of the free electrons. 


Due to the translational symmetry of the Hamiltonian of a two-electron ion, the degrees of free¬ 
dom reduce from nine to six by separating the motion of the centre of mass. These six coordinates 
can be taken as the sides of the triangle ri, r 2 , ri 2 formed by the three particles, he., two electrons 
and the hxed nucleus and the Eulerian angles (0, 0, 'ip) dehning the orientation of this triangle in 
space. The wave function obeying symmetry properties under particle exchange may be written 
as [S] 

« = 5] [/r (>- 1 , >- 2 , 9i 2) B2+(9, V>) + /r (>- 1 . >' 2 ,912) B2-(9, i>)] (6) 
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Qyi is the angle between rt and rf. The summation in eq. ( 6 ) goes over every alternate value of k, 
where k = \k\. k is the angular momentum quantum number about the body hxed axis of rotation 
whose value satishes k ^ L, L being the total angular momentum quantum number. The symmetric 
top functions D'l'^ and ^^e eigenfunctions of angular momentum operator of the two 

electrons. The rotational invariance of the Hamiltonian makes it possible to express the variational 
equation of two electrons in the held of a hxed nucleus in terms of three independent variables ri, r 2 
and ri 2 (or ^ 12 ). The reduction of the Eulerian angles from the variational equation is an immediate 
consequence of the spherical symmetry of the held. 


For state 


The variational equation for states originating from two s-electrons (Isns conhguration) is given 
by, following ref. na. 


s 


/ 


1 ' 1 ' frJ-rf + riA f fd^s\ 

2 \driJ 2\dr2j V 2riri2 J \ dn J \dr12 J 


- rg + 

2r2ri2 




+ {V-E)'^l 


dVr 


Ti,r2,ri2 


0 


(7) 


where the upper limit of integration for ri and r 2 is R in contrast to inhnity for the free atomic case. 
The upper and lower limits of integration for ri 2 are (ri + r 2 ) and |ri — r 2 |, respectively. The volume 
element is expressed as 


dVri,r2,ri2 = rir2ri2dridr2dri2 ( 8 ) 

For S'-states, L = k = 0 and the wave function T s can be written as 

Ss = /o“bS = fs + Js (9) 

where ''" 2 , ^ 2 ) = /s(?" 2 , ri 2 ). The correlated wave function [H] is written as 

/s(ri,r2,ri2) = (i?-ri)(i?-r2)/(ri,r2,ri2) (10) 


with 


/(ri,r 2 ,ri 2 ) = e EEE Cimnr[r^r^2 (H) 

/>0 m >0 n >0 

(Ti and (72 are the nonlinear parameters taking care of the ehect of radial correlation in the wave 
function whereas the angular correlation ehect is incorporated through diherent powers of ri 2 . C’s 
are the linear variational parameters. The total number of parameters {N) in the basis set is dehned 
as the total number of diherent sets (eq. 11) taken in the expansion of /(^i, r 2 , 712 ). 


For ^P° state 

For P state of odd parity (L = 1 and k = +1, —1), the total wave function can be written as [IB] 

Tp = fl+D\+ + (12) 

The origin of ^P° state due to sp conhguration of two-electron atoms can be shown by expressing 
DI~^ in terms of individual polar coordinates {9i, 0 i; 6 * 2 , ^ 2 ) as given below: 


0 , - 0 ) cos 


6^12 

'Y 


Pg (cos 6 'i)P°(cos ^ 2 ) + interchange 


(13) 
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Similar expressions for D\ can be derived. After integration over the Eulerian angles, the variational 
eqnation for ^P° states originating from Isnp conhguration reduces to [IH], 



1 

sin 6*12 



+ 2(1/-B){(/;+)V (/;-)'} 


dVr 


ri,r2,di2 


= 0 


The correlated functions /i’*' and fl can be written as 
= ifp + fp) cos ^ 


(14) 


(15) 


fl = ifp - fp) sin ^ (16) 

where 

/p(ri,r2,ri2) = (R - ri){R - r2)giri,r2,ri2) (17) 

and 

fi'(ri,r2,ri2) = EEE Dimnr{r^r^2 (18) 

/>1 m >0 n >0 

The nonlinear parameters i.e. cr’s in eq. (11) and p’s in eq. (18) for S and P states respectively 
are optimized separately using Nelder-Mead algorithm mi- The linear variational parameters i.e. 
CimnS and DimnS aloug with the energy eigenvalues are obtained by solving the generalized eigenvalue 
equation 

KC = EgC (19) 

where is the Hamiltonian matrix, ^ is the overlap matrix, C_ is the column matrix consisting of 
linear variational parameters and E is the corresponding energy eigenvalue. The wave function is 
normalized for each conhning radius R to account for the reorientation of charge distribution within 
the Wigner-Seitz sphere. All computations are carried out in quadruple precision to ensure better 
numerical stability for extended Hylleraas basis sets within a finite domain. 


The variational equation for nl {^L) states of one electron ion within ion-sphere of radius R can 
be written as 
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dl 

dr 


+ 


+ 1) s2 


r ) + {Visir) - E} 


r^dr = 0 


( 20 ) 


where, the one-particle effective potential Vjsir) is taken from eq. (2) with A'^g = 1. The radial 
function /(r) is given by 


/(r) = {R- ry 

i 


( 21 ) 
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where k = Q and 1 for ‘^S and states respectively. In this calculation, we have taken 21 dif¬ 
ferent nonlinear parameters (pi’s) in a geometrical sequence pi = Pi-i'y, 7 being the geometrical 
ratio [ISIIIH]. Such choice of non-linear parameters enable us to cover the full region of space in a 
flexible manner by adjusting 7. The energy values £”s and linear co-efficients Q’s are determined 
from eq. (19). 


The truncation of wavefunction at a hnite distance (eq. 3) imposes a thermodynamic pressure 
upon the ions which increases with increase of Ue inside the sphere. We have calculated the pressure 
felt by all the hydrogen-like and helium-like ions in their respective ground state using the hrst law 
of thermodynamics. However,for excited states having hnite lifetime, this approach is not valid as 
the equilibrium criteria is not maintained. Under an adiabatic approximation, the pressure on the 
ions in the ground state can be expressed as [TT] 


1 dE 


( 22 ) 


Evaluation of two-electron integrals 

The correlated two-electron basis integrals arising in the present calculations are of the form 






H(a, b, c; a, (3; R) = 




rl 


rri+r2 


rU dridr2dri2 


'In->'21 


r-R 


a „-ari I b „-/ 3 r 2 


r‘)'i+r2 


-R 


fr’2 


e 


To e 


rfn dridr2dri2 + r^e 


b „-pr 2 / „a „-ari 


rri+r2 


e 


'0 


'0 


ri2 dridr2dri2 

(23) 


For S states, a>0,6>0,c>0 while for higher angular momentum states {P, D etc.), integrals 
with a = — 1 also arises. After integration, the ri 2 part of eq. (23) can be expanded as 


^ [(ri + ri)”+‘ - (n - ri)”+‘] = ^ 


n -I- 1 


.^(2i+l)!(n-2»)! 


r” 2 ^^ 2 ^+l even] 


(24) 


For odd ‘n’, the upper limit of the sum in the right hand side would be replaced by The integrals 
from eq. (23) then reduce to the form 

ry 


e-^^dx = / e-^^dx - / e'^^dx = 


k\ 






i=o 


r- 


(25) 


A is a positive real number and A; is a non-negative integer and we have used the standard integral 


r-oo 7 I 

..k-Xx- bv. 


x^e-^^dx = 


yfc+1 


(26) 


The integral A(a, b, c; a, /3; R) is now evaluated for two different cases. 


Case I: a > 0, 6 > 0, c > 0 

An exact analytical expression for A{a,b, c] a, (3] R) corresponding to a > 0, 6 > 0, c > 0 has been 
derived in a straightforward way using eq. (25) and the numerical values are displayed in table-1. 
In the hrst column of table-1, different powers of ri, r 2 and ri 2 i.e. a, b and c are given. For each 
set of (a, b, c), the non-linear parameters (a, (3) given in the second column of table-1 are varied from 
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very low to high values as obtained from the optimized values corresponding to different cases in the 
present work. The values of R varied in a wide range for each set of (a, b, c) and (a, (3) are given 
in the fourth column of table-1. All the valnes of integrals are given in the in the last colnmn of 
table-1. We have checked the results with standard mathematical software {e.g. Maple) to ensnre 
the nnmerical stability of the expression for A(a, b, c; a, /3; R) over the complete range of R . 


Case II: a = —1, 6 > 0, c > 0 

After fnll expansion of the integral A(—1, b, c; a, /5; R) over ri 2 and r 2 by using eqs. (24) and (25) an 
integral I{a, /5; R) arise which takes the form 

rR Q-ari _ g-(a+/3)ri 

I(a,/3;R)= / - dri (27) 

Jo 

The above integral I{a, (3] R) is actnally a converging inhnite series with oscillatory terms. We have 
tested the evaluation of the term I{a,/3; R) in two different approaches. 


i) We can expand the exponential fnnctions to evalnate the integral as 



g-ari 


_ -{a+/3)ri 

— - dn 

n 


E 

q=0 

oo 

E 

q=l 


/o n 


-1)^ 


q.q\ 




{a® - (a -f /3)'^}r^ 


(a + /3)‘i 


dri 


(28) 


The expression (28) gives accnrate valne of integrals where the upper limit R is small, but fails to 

produce results when R is snfhciently high. 

ii) Alternatively, the integral I{a,l3] R) may be written as 



g-ari 


_ g-(a+/3)n 

—- dn = 

n 



(1 


— e dri = 

q=l 


(-l)<?-i/3<? 

q\ 


(‘R 


rf °‘^^dri 


(29) 


The ri-integral in the r.h.s of eqnation (29) is then evalnated using expression (25). 

The integral I{a,l3]R) is calculated by using both the expressions given in eq. (28) and eq. (29). 
All the results corresponding to different sets of {a, (3; R) are given in table 2 which shows excellent 
agreement among the resnlts except for some high values of R used in eq. (28). On the other hand, 
eq. (29) yields excellent results irrespective of the values of the parameters (a, /5) over the complete 
range of R. In eq. (28), a term R'^ appears in the nnmerator that increases with increase in q. 
For low valnes of R, this term is balanced by q\ in the denominator bnt for high R, a nnmerical 
instability appears because within the first few terms, R'^ bonnces more rapidly than q\. In contrast, 
a term appears in eq. (29) [after expanding the ri-integral according to eq. (25)] which falls 

rapidly as q increases dne to the presence of the exponential term. To have a better nnderstanding of 
the integrals, we have also checked the convergence of I{a,/3] R) evalnated using eqs. (28) and (29) 
by increasing the nnmber of terms in the inhnite series and displayed the convergence behavionr in 
table 3 for i? = 100.0 and 0.2 and two sets of {a,/3). It appears from table 3 that for R = 100.0 the 
valnes derived from eq. (28) are clearly not acceptable bnt for low R, the hnal results match exactly 
though the convergence is slow for eqnation (28). We have hnally used eq. (29) to calculate the 
energy eigenvalnes in the present work and taken 1000 terms in the corresponding inhnite series to 
ensnre the desired level of accnracy. In table 4 we have given the valnes of integral A(—1, b, c; a, /5; R) 
corresponding to diherent sets of parameters. We have further observed that the integrals [eq. (29)] 
corresponding to R = 100 yield same result as obtained by nsing eq. (26) for R = oo. This is mention 
fnrther that all the integrals are checked with standard mathematical software. 



3 Results and discussions 


The energy eigenvalues of He-like C, Al and Ar in Isns [n = 1 — 3] and Isn'p {^P°) [n' = 2 — 4] 
states have been calculated within SCP environment using the IS potential. We have studied the 
convergence of the energy values w.r.t. the number of terms {N) in the wave function. Table 5 
shows the convergence behaviour of in (^5'^) state for some selected values of R. We have 

obtained a similar convergence pattern for all the other ions and also for the excited states under 
consideration. The size of the basis has been extended systematically to = 161 and 149 for ^5'^^ 
and ^P° states respectively with / + m + n = 10 [eqs. (11) and (18)]. The convergence of the energy 
values are obtained at least up to the sixth signihcant digits. In fact, for some cases e.g. Is^ (^<5'^) 
state of with R = 0.47 a.u., we have obtained convergence of energy values up to the eighth 
decimal place, as is evident from table 5. The above observation ensures that the present method 
can deal with extended basis sets to yield sufficiently accurate energy values within a finite limit. 

The energy values of He-like C, Al and Ar in Isns (^S^) [n = 1 — 3] and Isn'p (^P°) [n' = 2 — 4] 
states within ion-sphere of different radii (i?) are displayed in tables 6, 7 and 8 respectively. We have 
also listed the energies of respective H-like ions in ns (^S) [n = 1 — 2] and n'p {^P) [n' = 2 — 3] 
states. It is worthwhile to mention that under one-component plasma approximation, the IS radius 
for a two-electron ion would differ from that for a one-electron ion corresponding to the same plasma 
electron density. We see that as Ue increases, the energy levels move towards continuum which is 
a clear manifestation of the positive nature of IS potential. To check the overall behavior of the 
results, we have plotted the energy values {—E) of bound Isns (^S*®) [n = 1 — 3] and Isn'p (^P°) 
[n' = 2 — 4] states of with respect to IS radius {R) in hgure 1. It is evident from hgure 1 
that the energy values remain almost unaltered for a range of R and after that rapidly approaches 
towards the destabilization limit. Hence the variation produces a PneP around some particular value 
of R. For higher excited states, this PneP appears at a higher value of R. All other ions also show 
same features. Similar behaviour of energy values of He-like ions inside a spherical impenetrable box 
(referred to as ‘Coulombic sphere’ hereinafter) was reported in a recent publication [H] where the 
potential inside the box was purely Coulombic. Within the ion-sphere, energy value of the positively 
charged ion is modihed for two factors: 

1. The environment envisaged by IS potential which is governed by plasma electrn density and 

2. The truncation of wave function at a hnite distance that generates a pressure on the system. 

In order to asses the effect each factor on the energy eigenvalues, we have also studied separately 
the modihcation of energy values of two-electron ions due to the truncation of the wave function 
at different radii of Coulombic sphere. The ground state energy of a ‘/rce’ ion where the wave 
function is inhnitely extended is —32.406247 a.u. whereas within a Coulombic sphere and ion-sphere 
both having a radius of 20.0 a.u, the energy values are —32.406247 and —31.806294 a.u. respectively. 
It shows that for a large box radius, almost 100% of the shift in the energy is due to the effect of 
plasma. The truncation of wave function becomes significant when the size of the sphere is reduced. 
At a radius of 0.7 a.u., the ground state energy values of ion within Coulombic sphere and 
ion-sphere are —31.192275 and —14.880628 a.u. respectively which shows an effect of almost 7% 
on the shift of energy level is coming from the truncation of wave function. This effect increases to 
22.4% and 26.2% for truncation radius of 0.5 and 0.4692 a.u. respectively. 

A closer look at the results quoted in tables 6 — 8 leads ns further to the following observations. 

1. Decrease in number of excited states: For two-electron ions Al^^^ and Ar^^^ we 

see that as rig increases, the ions become less bound and also the number of excited states 
decreases. For example, exists in the ground state up to i? = 0.4692 a.u. but ls2s (^(F®) 
state ceases to exist after R = 0.9017 a.u. and ls3s (^5''’) destabilizes after R = 1.3761 a.u. 
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Similar feature is observed for all the ions and also for ^P° states. For H-like ions of C, Al and 
ylr, the 2s state destabilizes much before Is with increase of R. 

2. Reduction of ionization potential: Ionization potential for a two-electron ion is dehned as 

the amount of energy required to ionize one electron from the ground state (Is^). It is observed 
from tables 6 — 8 for all the ions that with increase in plasma density, IP decreases and beyond 
certain density, the two-electron energy levels move above the one-electron threshold. We have 
studied the variation of IPD of two-electron ions w.r.t. Ue from the difference of IP within 
and without {i.e. free case) the surrounding plasma environment. In £gure-2, we have plotted 
the IP and IPD for Al^^+ as a function of plasma electron density. The energy required to 
ionize the outer electron from ls3p{^P°) state i.e. IP for ls3p{^P°) state of and the 

corresponding IPD w.r.t. Ue are also included in figure-2. The effect of surrounding plasma 
on different two-electron energy levels should be different and consequently, IPDs should differ 
from one conhguration to another. It is evident from £gure-2 that the present observation 
corroborates this fact. The two-electron states would give rise to spectral lines via dipole 
transition until they merge into the one-electron continuum. For example, table-7 shows that 
ls3p{^P°) state of can survive up to the density of 8.11 x 10^^/c.c. whereas it crosses 

the corresponding Is threshold after the density of 2.11 x 10^^/c.c and consequently, the Hep 
line originating from ls3p{^P°) —)■ ls^(^iS'®) emission is not expected to be observed after this 
density. In table 9, we have listed the critical electron densities after which different spectral 
lines of H-like and He-like Al disappear. The densities are calculated from IS radius according 
to both SP model and EK model of determining IPDs following eq. (4) and (5) respectively. 
For Lyp line of and Hep line of Al^^^, the present electron densities calculated by using 

SP model are in good agreement with experimental observation [33]. For disappearance of 

line, the only theoretical calculation of plasma density is due to Preston et al. [31] where 
a possible range of densities is given. No experimental result is available for comparison in 
this context. Our results obtained by using SP model of IPD indicate that the He.y line of 
Al^^^ would disappear after a plasma density of 5.0 x 10^^/c.c, as is given in table-9. We 
mention that the disappearance of both Lyp and Hep lines are experimentally observed at the 
density of 2.2 x 10^^/c.c. whereas the Lyp line should survive more than the Hep line. Present 
results along with ref. [35] as depicted in table-9 establish the fact explicitly. A more accurate 
experimental measurement is therefore necessary for proper plasma diagnostics. 

In an earlier experiment, Nantel et al. [22] observed Hca, Hep and He.y lines of at plasma 
density 1.5 x 10^^/c.c. In this experiment the densities corresponding to disappearance of such 
He-Iike lines are not explored. However, table-9 shows that He-like lines of (7^^ vanish well 
above the density of 1.5 x 10^^/c.c. Hence, the existence of such He-like lines of at the 
density 1.5 x 10^^/c.c. as observed by Nantel et al. [22] are consistent with present calculations. 
Similar comparisons have been done with other earlier experiments of Saemann et al. [23] and 
Woolsey et al. for spectral lines of Al^^~^ and Ar^^~^ respectively and the present results are 
in agreement with the experiments. Accurate measurement like the Orion laser experiment [33] 
is necessary to confirm the present theoretical predictions for disappearance of spectral lines of 

and 

3. Quasi-bound states of two-electron ions: Quasi-bound states or continuum bound states 
may be found in continuous part of the spectra for electronic conhnement under different po¬ 
tentials [H] and has also been observed experimentally [SU] . These states have great structural 
similarity with the discrete energy levels. For a two-electron ion, the ground state and all singly 
excited energy levels, in general, lie below the first ionization threshold. Tables 6-8 show that 
for high values ‘i?’ {i.e. almost free case), this feature is maintained for all the ions but as R 
decreases, all singly excited states of two-electron ions become less bound more rapidly than 
the respective one-electron ion. For example, at i? = 20.0 a.u. the energy values of as 
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reported in table-6 lie below the Is threshold of . At R = 5.0 a.u. the ls4p {}P°) state 
moves above the Is threshold but lies below 2s threshold. Similarly, at i? = 2.0 a.u. Is3s 
and ls3p {^P°) states lie above the Is threshold and below 2s threshold. At i? < 1.7 a.u. 
the 2s level of destabilizes and we observe well-converged (up to 7th signihcant digits) 
energy level of ls4p {^P°) state of embedded in one-electron continuum. Similar feature 
is obtained for other ions also and is being reported for the hrst time in SCP. 

4. Incidental degeneracy and level crossing: For a free two-electron ion, the energy value of 
ls2s state is more negative than ls2p {^P°) state. Tables 6-8 establish this fact for high 
values of R corresponding to all the two-electron ions. For example, the ls2s level of 
lies below ls2p (^P°) level for IS radius up to 7? = 2.0 a.u. At i? = 1.5 a.u., the ls2s 
state moves above the ls2p {^P°) level. These results show that an ^incidental degeneracy^ |13] 
has taken place for ls2s and ls2p {^P°) states of at some value of R between 1.5 

and 2.0. a.u. and then a devel crossing^ occurs between two states having different symmetry 
properties. The phenomenon of incidental degeneracy was reported in case of shell-conhned 
hydrogen atom by Sen |13] where two initially non-degenerate states are brought to a same 
energy level by adjusting external parameters. For a two-electron ion, we report incidental 
degeneracy for the hrst time within SCP environment. After the level crossing, the ls2s 
state of destabilizes (at R = 0.9017 a.u.) much before ls2p CP°) state (at R = 0.797 a.u.). 
Similarly, for ls3s (^iS®) and ls3p i^P°) states of incidental degeneracy and subsequent 
level crossing are observed at a value of R lying somewhere between 5.0 and 3.0 a.u. We 
observe similar phenomena for other ions also. For 2s and 2p states of H-like ions embedded 
in SCP, incidental degeneracy' and devel crossing' phenomena are evident from tables 6 — 8 
and being reported here for the hrst time. 

In table 10, we have displayed the comparison of present results for transition energies of , Al^^^ 
and Ar^^~^ with experimental measurements [22H21] and existing theoretical calculations |3S]- Table- 
10 clearly depicts that the present results are in good agreement with the experimental measurements 
as compared to other theoretical results. The diherence between theoretical and experimental results 
is due to the relativistic ehects which is more signihcant for Al^^~^ and Ar^^^ than For 
conhned in ion-sphere, Sil et al. [36] have made a comparison of theoretically calculated transition 
energies with the experiments but the plasma electron density reported by them is one order less 
than the experimental condition of 1.0 x 10^^ per c.c. [23] and the IS radius of 2.057 a.u. [36] also 
does not correspond to any of the densities. Actually, according to eq. (4), IS radius of 2.057 a.u. 
corresponds to a plasma electron density of 2.96 x 10^^ per c.c. for We have calculated the 

energies corresponding to Is^ Isnp [n = 2 — 4] transitions of Ar^^~^ for both experimental [23] 
and theoretical [36] conditions and given the values in table-10. It is seen from tabIe-10 that the 
Is^ —?• Is5p transition energies estimated from present results are more away from experiments as 
compared to other transitions. We note that present size of the trial wave function for a high-lying 
state like ls5p i^P°) is not sufficient and needs to be increased to obtain greater accuracy. 

The thermodynamic pressure experienced by H-like and He-like C, Al and Ar in their respective 
ground states have been calculated for different values of IS radius R using eq. (22) and the results 
are given in tabIe-11. It is clear from table-11 that as n^ increases, the pressure upon the ion increases 
and the ion moves towards destabilization. We observe that for a low value of Ug, the pressure upon 
the one-electron ion is higher than the respective two-electron ion and after a certain increase of 
Ue, the pressure on the two-electron ion exceeds the pressure experienced by the corresponding one- 
electron ion. With a view to studying the variation of thermodynamic pressure (P) with respect to 
the IS volume (H) under an adiabatic expansion, we have tried to fit the results for the two-electron 
ions obtained from the present calculations according to the ideal gas relation, 

PW = constant or InP + 'ylnV = constant (30) 
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where 7 is the ratio of two specihc heats. From a least square ht of In P vs. In V plot, the value of 
7 comes out to be close to 1.4 for all the two-electron ions. To be precise, for and , 

the values of 7 are 1.41, 1.37 and 1.37 respectively. 

4 Conclusion 

Accurate analytical evaluation of the two-electron correlated integrals in Hyllerass coordinates within 
a hnite limit has been performed. The intricacies of such calculations have been discussed in detail 
and the general applicability of these integrals has been established for arbitrary values of physi¬ 
cally acceptable parameters. This methodology has immense potential to be nsefnl for evalnation 
of the energy valnes and other spectral properties for three-body ionic and exotic systems placed 
within different external conhnements snch as strongly and moderately conpled plasma, fullerene 
cages, barrier potential, potential well etc. With the recent advancement in experimental techniqne, 
the present methodology becomes relevant for calculating accnrate plasma electron density from 
the spectral analysis of hydrogen and helinm-like ions. We conclnde that the ion-sphere potential 
where the electron density is calculated by using the SP model of IPDs provides a realistic pic- 
tnre of ions embedded in SCP environment. The present non-relativistic results reported here can 
be nsefnl for plasma diagnostics and the non-relativistic energy valnes can serve as benchmark for 
fnture calcnlations to estimate relativistic and QED effects on two-electron ions within hnite domain. 
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Ion sphere radius ( R) in a.u. 


Figure 1: Variation of eigen energies {—E) of bound [n = 1 —3] and lsn'p{^P°) [n' = 2 — 4] 

states of with respect to ion-sphere radius R. 



Figure 2: Variation of IP and IPD for (^5'^) and ls3p(^P°) states of with respect to plasma 
electron density. 
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Table 1: Values of integral A{a,h, c; a, /3; R) with a > 0,& > 0,c > 0. The notation x{y) indicates 

X X IQJ'. 


( a , b , c ) 

a 

P 

R 

A { a , b , c ; a , / 3 ; R ) 


( 0 , 0 , 0 ) 

0.62450527 

0.41287135 

100.0 

0.7477 2634 8987 8847 

(+ 01 ) 




2.0 

0.1578 7453 6391 8980 

(+ 01 ) 




0.2 

0.4688 3009 2327 4260 

(- 02 ) 


8.92934001 

5.97270373 

100.0 

0.2516 4820 0437 6827 

(- 02 ) 




2.0 

0.2516 4546 4241 1431 

(- 02 ) 




0.2 

0.9770 3676 6565 7275 

(- 03 ) 


17.42010556 

10.32300145 

100.0 

0.4008 8345 0596 5748 

(- 03 ) 




2.0 

0.4008 8344 9906 8468 

(- 03 ) 




0.2 

0.3006 0461 4922 4955 

(- 03 ) 

( 2 , 3 , 1 ) 

0.62450527 

0.41287135 

100.0 

0.1578 2446 3158 5587 

(+ 06 ) 




2.0 

0.9755 4412 5697 4906 

(+ 01 ) 




0.2 

0.4326 2415 0573 9082 

(- 07 ) 


8.92934001 

5.97270373 

100.0 

0.5960 1910 9088 6659 

(- 05 ) 




2.0 

0.5913 0152 3072 3085 

(- 05 ) 




0.2 

0.4824 7810 4677 3954 

(- 08 ) 


17.42010556 

10.32300145 

100.0 

0.2667 8377 6884 8479 

(- 07 ) 




2.0 

0.2667 8110 4519 4153 

(- 07 ) 




0.2 

0.7198 6505 5902 0197 

(- 09 ) 

( 3 , 4 , 6 ) 

0.62450527 

0.41287135 

100.0 

0.5959 0474 3300 4562 

(+ 14 ) 




2.0 

0.3133 2861 0504 5369 

(+ 04 ) 




0.2 

0.1612 2935 4296 9925 

(- 11 ) 


8.92934001 

5.97270373 

100.0 

0.1741 4651 3311 5703 

(- 04 ) 




2.0 

0.1311 5563 4874 3337 

(- 04 ) 




0.2 

0.1427 2637 0766 0009 

(- 12 ) 


17.42010556 

10.32300145 

100.0 

0.1312 3966 9657 3856 

(- 08 ) 




2.0 

0.1305 9386 1401 5103 

(- 08 ) 




0.2 

0.1601 5416 4752 0373 

(- 13 ) 


Table 2: Values of integral I{a,l3]R). Results obtained by using eq. (28) and eq. (29) are given in 
consecutive rows respectively. The notation x{y) indicates x x 10^. 


a 

P 

R 

I { a , p ; R ) 

0.62450527 

0.41287135 

100.0 

0.2910 1096 5167 7626 (+ 9 ) 
0.5074 9055 6270 2974 



2.0 

0.4049 6455 3461 7721 

0.4049 6455 3461 7721 



0.2 

0.0760 8455 8526 4357 

0.0760 8455 8526 4357 

8.92934001 

5.97270373 

100.0 

0.2617 9165 1879 3351 (+ 602 ) 
0.5121 5588 2201 4390 



2.0 

0.5121 5588 1268 8812 

0.5121 5588 1268 8812 



0.2 

0.4595 0838 7239 3843 

0.4595 0838 7239 3843 

17.42010556 

10.32300145 

100.0 

0.9419 3093 2537 1151 (+ 872 ) 
0.4653 6238 2193 5131 



2.0 

0.4653 6238 2193 5052 

0.4653 6238 2193 5131 



0.2 

0.4588 5876 8812 0597 

0.4588 5876 8812 0597 
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Table 3: Convergence of the integral I{a,l3-,R) w.r.t. the number of terms {q) in the inhnite series 
using equations (28) and (29). The notation x{y) indicates x x 10^. 


a 

/3 

R 

Q 

I { a , j 3 ] R ) 





Eq . ( 28 ) 



Eq . ( 29 ) 


0.62450527 

0.41287135 

100 

10 

- 0.3568 5265 8756 6828 

(+ 13 ) 

0.5068 

9300 

3143 

0272 




20 

- 0.3554 0371 2221 5343 

(+ 21 ) 

0.5074 

8564 

3805 

4664 




50 

- 0.2755 0482 1577 0421 

(+ 35 ) 

0.5074 

9055 

6262 

6217 




100 

- 0.2124 7555 4344 8736 

(+ 42 ) 

0.5074 

9055 

6270 

2974 




1000 

0.2910 1096 5167 7626 

(+ 09 ) 

0.5074 

9055 

6270 

2974 



0.2 

10 

0.0760 8455 8526 4356 


0.0760 

8455 

8526 

4357 




20 

0.0760 8455 8526 4357 


0.0760 

8455 

8526 

4357 




50 

0.0760 8455 8526 4357 


0.0760 

8455 

8526 

4357 




100 

0.0760 8455 8526 4357 


0.0760 

8455 

8526 

4357 




1000 

0.0760 8455 8526 4357 


0.0760 

8455 

8526 

4357 

17.42010556 

10.32300145 

100 

10 

- 0.7343 5142 7504 1944 

(+ 27 ) 

0.4651 

7579 

3447 

4452 




20 

- 0.1488 0713 2757 9564 

(+ 50 ) 

0.4653 

6186 

8073 

0945 




50 

- 0.9285 2359 3991 9110 

(+ 106 ) 

0.4653 

6238 

2193 

4812 




100 

- 0.2137 8120 9553 6411 

(+ 185 ) 

0.4653 

6238 

2193 

5131 




1000 

- 0.2609 9042 7983 5017 

(+ 873 ) 

0.4653 

6238 

2193 

5131 



0.2 

10 

0.2166 5736 2229 1603 


0.4588 

5852 

3846 

0869 




20 

0.4588 5558 5227 7264 


0.4588 

5876 

8812 

0596 




50 

0.4588 5876 8812 0597 


0.4588 

5876 

8812 

0597 




100 

0.4588 5876 8812 0597 


0.4588 

5876 

8812 

0597 




1000 

0.4588 5876 8812 0597 


0.4588 

5876 

8812 

0597 
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Table 4: Values of integral A{a, b, c; a, 13] R) with a = —1, b > 0, c > 0. The notation x{y) indicates 

X X 10^. 


(a, b, c) a 


P R A{—l,b,c-,a,p-,R) 


(-1,0,0) 0.62450527 

8.92934001 

17.42010556 

(-1,3,0) 0.62450527 

8.92934001 

17.42010556 

(-1,6,4) 0.62450527 

8.92934001 

17.42010556 


0.41287135 100.0 

2.0 
0.2 

5.97270373 100.0 

2.0 
0.2 

10.32300145 100.0 

2.0 
0.2 

0.41287135 100.0 

2.0 
0.2 

5.97270373 100.0 

2.0 
0.2 

10.32300145 100.0 

2.0 
0.2 

0.41287135 100.0 

2.0 
0.2 

5.97270373 100.0 

2.0 
0.2 

10.32300145 100.0 

2.0 
0.2 


0.5954 2692 0365 7702 (-H01) 
0.2328 7535 9564 1012 (-hOl) 
0.5410 3409 5795 4025 (-01) 
0.2871 3769 8867 2015 (-01) 
0.2871 3526 4943 5828 (-01) 
0.1630 9271 2387 9263 (-01) 
0.8733 9216 7581 8471 (-02) 
0.8733 9216 6380 3946 (-02) 
0.7244 1513 1903 8541 (-02) 
0.6567 2317 0086 4201 (-e03) 
0.4615 4100 9066 9881 (-H01) 
0.1352 9814 0488 3989 (-03) 
0.1048 5277 9488 0441 (-02) 
0.1046 0034 4655 7635 (-02) 
0.2845 3033 4471 0229 (-04) 
0.6034 8938 8061 1726 (-04) 
0.6034 8827 2996 3596 (-04) 
0.8940 2524 7131 8926 (-05) 
0.1994 0115 8159 1710 (-^12) 
0.3251 2782 4902 2441 (-h03) 
0.1201 4454 0849 6434 (-08) 
0.2402 8418 0970 4088 (-02) 
0.1563 9297 5468 9225 (-02) 
0.1788 1043 5252 0599 (-09) 
0.2982 9940 1992 8102 (-05) 
0.2960 5934 0875 9482 (-05) 
0.4140 7243 9026 4580 (-10) 


Table 5: Convergence of energy values {—E a.u.) of state of with respect to number of 

terms {N) in wave function within Ion-sphere radius R a.u. 


State 

N 


—E for two-electron 

ions 


77 = 20.0 

0.7 

0.5 

0.47 

0.4692 

ls2(i^^) 

13 

31.8060 7622 

14.8798 5637 

3.4091 6161 

0.1035 1773 

0.0053 4000 


22 

31.8062 6559 

14.8805 9032 

3.4092 5286 

0.1037 0540 

0.0055 3914 


34 

31.8062 9082 

14.8806 1939 

3.4092 6383 

0.1037 1810 

0.0055 5223 


50 

31.8062 9351 

14.8806 2562 

3.4092 6610 

0.1037 1980 

0.0055 5423 


70 

31.8062 9412 

14.8806 2707 

3.4092 6664 

0.1037 2017 

0.0055 5462 


95 

31.8062 9431 

14.8806 2746 

3.4092 6678 

0.1037 2026 

0.0055 5471 


125 

31.8062 9439 

14.8806 2759 

3.4092 6682 

0.1037 2029 

0.0055 5474 


161 

31.8062 9443 

14.8806 2763 

3.4092 6683 

0.1037 2029 

0.0055 5475 
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Table 6: Energy eigenvalues {—E a.u.) of (7^^ and within ion-sphere of Radius R a.u. Densities (per c.c.) are determined from radius by 
using eq. (4). The uncertainty of the calculated energy values is of the order of 10“® a.u. 


Plasma 

density 



-E 

for two-electron ions 




—E for one-electron ions 


R 

1^2 

ls2s 

ls3s 

ls2p 

ls3p 

Isip 

R 

Els 

E2s 

E2p 

Esp 


CA+ 







C«+ 





8.05(20) 

20.0 

31.806 294 

20.622 433 

18.819 815 

20.493 660 

18.782 365 

18.183 990 

21.544 

17.651 901 

4.152 167 

4.151 006 

1.298 528 

1.50(21) 

16.256 

31.668 147 



20.355 753 

18.645 771 

18.051 134 

17.511 

17.571 741 

4.072 241 

4.070 983 

1.216 051 

6.44(21) 

10.0 

31.206 629 

20.025 353 

18.234 022 

19.895 956 

18.195 421 

17.626 905 

10.772 

17.303 930 

3.806 085 

3.804 210 

0.941 359 

1.88(22) 

7.0 





17.710 079 

17.127 014 

7.540 

17.005 859 

3.512 120 

3.508 790 

0.638 130 

5.15(22) 

5.0 

30.009 311 

18.848 891 

17.115 098 

18.714 486 

17.073 177 

16.164 573 

5.386 

16.608 855 

3.126 251 

3.119 111 

0.239 771 

8.32(22) 








4.591 



2.885 861 

0.001 097 

1.01(23) 

4.0 






14.979 005 

4.309 

16.262 009 

2.796 224 

2.783 284 


2.39(23) 

3.0 

28.420 442 

17.332 928 

15.017 054 

17.183 256 

15.108 681 

12.307 491 

3.232 

15.685 381 

2.266 340 

2.237 178 


4.12(23) 

2.5 



13.431 934 


13.720 854 

9.562 152 

2.693 

15.225 721 

1.855 508 

1.812 626 


6.05(23) 

2.2 






6.922 023 

2.370 

14.850 960 

1.512 931 

1.468 923 


8.05(23) 

2.0 


15.256 681 

10.320 034 

15.180 308 

11.114 072 

4.454 707 

2.154 

14.539 678 

1.207 836 

1.175 386 


1.10(24) 

1.8 



8.209 279 



1.315 545 

1.939 

14.160 665 

0.797 941 

0.799 140 


1.31(24) 

1.7 



6.840 488 


8.309 520 

0.736 221 

1.831 

13.938 515 

0.529 722 

0.563 106 


1.57(24) 

1.6 

24.999 764 




7.004 681 


1.723 

13.689 353 

0.200 592 

0.283 113 


1.58(24) 

1.596 






0.030 797 

1.719 

13.678 753 

0.186 970 

0.271 734 


1.85(24) 








1.632 



0.003 072 


1.91(24) 

1.5 


12.478 314 

3.173 818 

12.790 730 



1.616 

13.407 990 




2.35(24) 

1.4 



0.691 166 


3.506 751 


1.508 

13.087 774 




2.38(24) 

1.3937 



0.516 288 




1.501 

13.066 138 




2.47(24) 

1.3761 



0.015 255 




1.482 

13.004 619 




2.93(24) 

1.3 


10.292 255 



1.124 349 


1.400 

12.720 153 




3.20(24) 

1.263 





0.096 911 


1.360 

12.569 908 




3.73(24) 

1.2 


8.722 057 


9.902 622 



1.293 

12.293 555 




6.44(24) 

1.0 

20.674 940 

3.865 288 


6.473 522 



1.077 

11.190 638 




7.52(24) 

0.95 


2.091 381 





1.023 

10.841 350 




8.55(24) 

0.91 


0.434 054 





0.980 

10.532 355 




8.79(24) 

0.9017 


0.059 416 





0.971 

10.464 484 




8.84(24) 

0.9 




3.864 058 



0.969 

10.450 413 




1.25(25) 

0.8 




0.209 360 



0.862 

9.499 773 




1.27(25) 

0.797 




0.077 557 



0.858 

9.466 811 




1.88(25) 

0.7 

14.880 628 






0.754 

8.198 121 




5.15(25) 

0.5 

3.409 267 






0.539 

2.956 871 




6.21(25) 

0.47 

0.103 720 






0.5063 

1.497 759 




6.23(25) 

0.4695 

0.042 434 






0.5057 

1.470 782 




6.24(25) 

0.4692 

0.005 555 






0.5054 

1.454 571 






Table 7: Energy eigenvalues {—E a.u.) of and within ion-sphere of Radius R a.u. Densities (per c.c.) are determined from radius 

by using eq. (4). The uncertainty of the calculated energy values is of the order of 10“® a.u. 


Plasma 

density 



—E for two-electron 

ions 




-E 

for one-electron ions 


R 

ls2 

1^25 

ls3s 

ls2p 

ls3p 

IsAp 

R 

Els 

E2s 

E2p 

Esp 


aA+ 












2.21(21) 

20.0 

159.382 029 

101.075 286 

90.920 760 

100.751 849 

90.826 863 

87.343 470 

20.589 

83.625 742 

20.250 901 

20.250 869 

8.508 494 

1.77(22) 

10.0 

157.732 211 

99.426 741 

89.277 651 

99.102 949 

89.183 005 

85.713 789 

10.294 

82.751 557 

19.377 826 

19.377 385 

7.638 134 

5.17(22) 

7.0 

156.318 323 

98.015 642 

87.878 476 




7.206 

82.002 366 

18.631 067 

18.629 929 

6.898 570 

1.00(23) 

5.617 

155.158 153 



96.533 764 

86.640 729 

83.249 437 

5.782 

81.387 603 

18.019 768 

18.017 402 

6.297 680 

1.42(23) 

5.0 




95.811 753 

85.932 485 

82.582 105 

5.147 

81.003 693 

17.638 853 

17.635 633 

5.926 068 

5.00(23) 

3.285 

150.992 205 



92.391 567 

82.632 916 

79.441 867 

3.382 

79.179 940 

15.840 819 

15.830 292 

4.207 565 

6.56(23) 

3.0 

150.039 702 

91.786 850 

81.867 981 




3.088 

78.675 129 

15.347 266 

15.332 246 

3.746 670 

1.00(24) 

2.607 

148.385 508 



89.814 745 

80.216 318 


2.684 

77.798 403 

14.495 335 

14.473 891 

2.966 481 

2.21(24) 

2.0 

144.557 994 

86.435 075 

76.751 134 

86.060 502 

76.632 056 

70.991 150 

2.058 

75.769 505 

12.554 606 

12.501 889 

1.209 389 

3.31(24) 








1.8 



11.309 886 

0.084 801 

5.25(24) 

1.5 

139.093 636 

81.216 562 

70.476 891 

80.776 714 

70.672 695 

59.958 894 

1.544 

72.871 972 

9.877 422 

9.760 672 


1.77(25) 

1.0 

128.240 365 

70.674 599 

51.222 379 

70.269 488 

53.774 530 

27.642 101 

1.029 

67.112 431 

4.647 358 

4.445 998 


2.43(25) 

0.9 




66.602 776 


13.910 306 

0.926 

65.206 039 

2.699 333 

2.596 821 


3.10(25) 

0.83 



34.498 851 

63.390 834 


1.358 983 

0.854 

63.604 555 

0.880 652 

0.952 002 


3.21(25) 

0.82 






0.585 110 

0.844 

63.354 014 

0.578 268 

0.692 836 


3.31(25) 

0.812 






0.088 777 

0.836 

63.149 234 

0.327 244 

0.479 227 


3.46(25) 

0.8 


61.584 852 

30.304 265 




0.823 

62.834 616 


0.011 954 


5.17(25) 

0.7 

114.496 531 

54.024 877 

11.937 344 

55.190 310 

21.877 996 


0.721 

59.806 678 




5.64(25) 

0.68 



7.201 566 




0.700 

59.097 854 




6.16(25) 

0.66 



2.007 102 




0.6794 

58.347 614 




6.18(25) 

0.6594 



1.843 590 




0.6788 

58.324 454 




6.37(25) 

0.65283 



0.022 248 




0.672 

58.067 871 




6.45(25) 

0.65 





11.929 838 


0.669 

57.955 830 




7.81(25) 

0.61 





2.170 876 


0.628 

56.265 059 




8.08(25) 

0.603 





0.262 866 


0.621 

55.946 962 




8.10(25) 

0.6025 





0.124 062 


0.6202 

55.923 954 




8.11(25) 

0.6021 





0.012 689 


0.6198 

55.905 567 




1.41(26) 

0.5 

96.523 118 

22.189 554 


30.281 348 



0.515 

50.268 732 




2.08(26) 

0.44 


2.370 513 





0.453 

45.729 585 




2.14(26) 

0.436 


0.724 569 





0.449 

45.379 873 




2.15(26) 

0.4351 


0.347 428 





0.448 

45.300 161 




2.16(26) 

0.43434 


0.026 971 





0.447 

45.232 635 




2.77(26) 

0.4 




3.174 966 



0.412 

41.877 044 




2.87(26) 

0.395 




1.274 781 



0.407 

41.331 971 




2.95(26) 

0.3918 




0.020 739 



0.403 

40.974 430 




6.56(26) 

0.3 

49.123 982 






0.309 

26.244 361 




1.28(27) 

0.24 

8.502 491 






0.247 

6.577 159 




1.419(27) 

0.232 

0.362 136 






0.239 

2.673 160 




1.424(27) 

0.2317 

0.038 819 






0.238 

2.518 185 






Table 8: Energy eigenvalues {—E a.u.) of and within ion-sphere of Radius R a.u. Densities (per c.c.) are determined from radius 

by using eq. (4). The uncertainty of the calculated energy values is of the order of 10“® a.u. 


Plasma 

density 



— 

E for two-electron ions 




-E 

for one-electron ions 


R 

1^2 

Is2s 

ls3s 

ls2p 

ls3p 

Isip 

R 

Els 

E2s 

E2p 

Esp 


aA^+ 







aA^+ 





3.22(21) 

20.0 

310.507 205 

196.041 306 

175.755 299 

195.578 034 

175.620 951 

168.619 605 

20.408 

160.750 518 

39.250 638 

39.250 583 

16.751 038 

2.58(22) 

10.0 

308.107 341 

193.642 370 

173.360 319 

193.178 832 

173.225 396 

166.234 270 

10.204 

159.501 078 

38.002 041 

38.001 721 

15.505 423 

2.06(23) 

5.0 

303.308 425 

188.850 886 

168.600 563 

188.385 210 

168.461 043 

161.553 149 

5.102 

157.002 621 

35.510 327 

35.507 887 

13.037 578 

9.55(23) 

3.0 






155.654 881 

3.061 

153.672 785 

32.208 489 

32.196 858 

9.835 118 

2.11(24) 

2.3 






151.183 198 

2.347 

151.140 923 

29.720 258 



2.96(24) 

2.057 

289.590 037 



174.748 078 

155.284 420 

148.653 416 

2.099 

149.859 811 

28.470 821 

28.436 728 

6.367 203 

3.22(24) 

2.0 

288.926 557 

174.593 816 

154.885 258 

174.092 315 

154.670 904 

147.924 331 

2.041 

149.514 344 

28.135 153 

28.099 008 

6.066 598 

7.64(24) 

1.5 



147.515 430 


147.294 212 

137.296 178 

1.531 

145.362 038 

24.149 851 

24.065 320 

2.548 088 

1.22(25) 








1.309 



21.374 828 

0.027 313 

1.49(25) 

1.2 



139.209 711 



123.092 431 

1.224 

141.218 060 

20.283 441 

20.102 175 


2.58(25) 

1.0 

265.062 279 

151.634 784 

129.063 495 

150.887 572 

129.848 525 

105.387 057 

1.020 

137.084 181 

16.548 570 

16.250 296 


5.03(25) 

0.8 



109.499 582 



72.434 154 

0.816 

130.907 580 

11.023 329 

10.593 320 


7.51(25) 

0.7 

244.788 710 

140.318 457 

92.135 024 

131.230 710 

97.682 171 

44.281 500 

0.714 

126.516 643 

6.834 684 

6.481 745 


9.39(25) 

0.65 






25.066 379 

0.663 

123.824 592 

4.030 093 

3.846 749 


1.13(26) 

0.61 






6.173 887 

0.622 

121.360 487 

1.249 801 

1.311 525 


1.18(26) 

0.6014 






1.639 257 

0.614 

120.788 926 

0.571 762 

0.762 695 


1.19(26) 

0.6 



64.536 211 




0.612 

120.694 413 

0.458 366 

0.622 378 


1.201(26) 

0.5987 






0.228 118 

0.6109 

120.606 199 

0.352 195 

0.544 660 


1.203(26) 

0.59843 






0.067 636 

0.6106 

120.587 871 

0.330 097 

0.523 396 


1.24(26) 








0.604 



0.048 119 


2.06(26) 

0.5 

218.150 195 

99.801 076 

17.342 729 

102.188 104 

37.003 115 


0.510 

112.615 414 




2.19(26) 

0.49 



10.878 182 




0.500 

111.632 327 




2.33(26) 

0.48 



3.985 508 




0.4898 

110.609 857 




2.344(26) 

0.4791 



3.342 920 




0.4889 

110.515 854 




2.345(26) 

0.47902 



3.285 614 




0.4888 

110.507 440 




2.41(26) 

0.47455 



0.035 558 




0.4842 

110.035 148 




2.83(26) 

0.45 





8.128 548 


0.4592 

107.280 138 




3.053(26) 

0.4387 





0.218 634 


0.448 

105.912 248 




3.059(26) 

0.43843 





0.022 169 


0.447 

105.878 796 




4.03(26) 

0.4 


63.058 578 


72.434 796 



0.408 

100.662 918 




6.01(26) 

0.35 


30.887 328 





0.357 

92.198 308 




7.87(26) 

0.32 


3.274 140 





0.326 

85.821 161 




8.01(26) 

0.318 


1.126 899 





0.3245 

85.350 316 




8.04(26) 

0.3176 


0.692 178 





0.3241 

85.255 417 




8.05(26) 

0.31751 


0.594 121 





0.3240 

85.233 983 




8.09(26) 

0.31698 


0.014 841 





0.323 

85.107 902 




9.55(26) 

0.3 

156.809 245 






0.306 

80.799 235 




1.08(27) 

0.288 




0.703 578 



0.294 

77.402 624 




1.09(27) 

0.2874 




0.09 253 



0.293 

77.224 113 




3.22(27) 

0.2 

64.368 049 






0.204 

34.699 257 




4.42(27) 

0.18 

25.307 474 






0.184 

15.611 535 




5.15(27) 

0.171 

2.741 398 






0.1745 

4.631 691 




5.24(27) 

0.1701 

0.269 467 






0.1736 

3.429 550 




5.25(27) 

0.17001 

0.019 935 






0.1735 

3.308 080 






Table 9: Critical plasma electron densities after which spectral lines of hydrogen-like and helinm-like 
Al disappear. Densities are obtained from IS radii according to both SP model [eq. (4)] and EK 
model[eq. (5)] for IPDs. The notation x{y) indicates x x 10^. 




Critical plasma electron density (per c.c.) 

z 

Spectral 

Present resnlts 

Other resnlts 


line 

SP modeP 

EK modeh 

Experiment 

Theory 

6 

^Va 

1.85(24) 

1.19(24) 




Lyp 

8.32(22) 

2.77(22) 





8.05(23) 

1.61(23) 




Hey 

5.15(22) 

1.03(22) 




He^ 

1 .88(22) 

3.76(21) 



13 

Lya 

3.52(25) 

2.71(24) 




Lyy 

3.31(24) 

2.55(23) 

2.2(24)= 

2.93(24)'^ 

2.64(24) - 3.3(24)= 


Hea 

2.43(25) 

2.53(24) 




Hey 

2.21(24) 

1.84(23) 

2.2(24)= 

2.44(24)-^ 

1.98(24) - 2.64(24)= 


He^ 

5.00(23) 

4.16(22) 


6.60(22) - 1.32(23)= 

18 

Lya 

1.24(26) 

6.90(25) 




Lyy 

1.22(25) 

6.78(23) 




Hcq, 

7.51(25) 

4.42(25) 




Hey 

7.64(24) 

4.49(23) 




He^ 

2.11(24) 

1.25(23) 




a: Ref. [55]; b: Ref. [30]; c: Ref. [33]; d: Ref. [3S]; e: Ref. [3T] . 
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Table 10: Transition energies (in A) in strongly coupled plasma environment. Densities are calculated 
according eq. (4). Conversion factors : (Ai?)^ = 12395/(Ai?)ev and 1 a.u. of energy = 27.21138 
eV. The notation x{y) indicates x x 10^. 



Ion-sphere 

Plasma density 

Transition 

Present 

Experimental 

Other theory”' 

Ion 

radius (a.u) 

(per c.c.) 

scheme 

results 

results 

Non-rel. 

Rel. 


16.256 

1.5(21) 

Is^ —)■ ls2p 

40.266 

40.268“ 

40.208 

40.700 




—)■ ls3p 

34.979 

34.998“ 

34.953 

35.325 




—)■ ls4p 

33.451 

33.469“ 

33.431 

33.773 




—)■ ls5p 

32.820 

32.773“ 

32.800 

33.133 


2.607 

1.0(24) 

Is^ ^ ls2p 

7.777 

7.75^ 

7.778 

7.774 




ls3p 

6.682 


6.684 

6.682 


3.285 

5.0(23) 

Is^ ^ ls2p 

7.773 


7.773 

7.770 




—)■ ls3p 

6.663 


6.664 

6.661 




—)■ IsAp 

6.366 



6.385 


5.617 

1.0(23) 

Is^ —)■ ls2p 

7.770 


7.769 

7.767 




—)■ ls3p 

6.648 

6.63'’ 

6.648 

6.646 




—>• IsAp 

6.334 

6.31'’ 

6.334 

6.332 




ls5p 

6.211 

6.17'’ 

6.205 

6.212 

^^16+ 

3.856 

4.5(23) 

Is^ ^ ls2p 

3.379 

3.363“ 




3.119 

8.5(23) 

Is^ ^ ls2p 

3.3811 

3.364“ 




3.061 

9.0(23) 

Is^ —)■ ls2p 

3.3813 

3.365“ 




2.955 

1.0(24) 

Is^ —)■ ls2p 

3.964 

3.984“ 






—)■ ls3p 

3.382 

3.365“ 






—)■ ls4p 

3.225 

3.168“ 




2.057 

2.96(24) 

Is^ —)■ ls2p 

3.966 


3.964 

3.950 




—)■ ls3p 

3.392 


3.382 

3.370 




^ IsAp 

3.232 


3.226 

3.215 


a: Ref. [52]; b: Ref. [23]; c: Ref. [22]; d: Ref. [3S] 
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Table 11; Thermodynamic pressure on the ground state of one and two-electron ions within ion- 
sphere. Conversion factor : 1 a.u. of pressure = 2.9421912(13) Pa. The notation x{y) indicates 


X X 10^. 


Plasma 

density 

8.05(20) 

6.44(21) 

5.15(22) 

2.39(23) 

1.57(24) 

6.44(24) 

1.88(25) 

5.15(25) 

6.21(25) 

6.23(25) 

6.24(25) 


Pressure (Pa) Plasma Pressure (Pa) Plasma Pressure (Pa) 


C^+ _ C^+ I density 

0.1755(09) 0.8778(09) 2.21(21) 0.4829(09) 

0.2807(10) 0.7019(10) 1.77(22) 0.7725(10) 

0.4478(11) 0.5603(11) 5.17(22) 0.3217(11) 

0.3432(12) 0.2581(12) 2.21(24) 0.4806(13) 

0.4130(13) 0.1666(13) 1.77(25) 0.7580(14) 

0.2779(14) 0.6681(13) 5.17(25) 0.3095(15) 

0.1540(15) 0.1794(14) 1.41(26) 0.1178(16) 

0.9278(15) 0.1772(14) 6.56(26) 0.1228(17) 

0.1297(16) 0.9241(14) 1.28(27) 0.3912(17) 

0.1304(16) 0.1123(15) 1.419(27) 0.4678(17) 

0.1309(16) 0.1164(15) 1.424(27) 0.4710(17) 


density 

0.2414(10) 3.22(21) 0.7024(09) 0.3512(10) 

0.1931(11) 2.58(22) 0.1123(11) 0.2809(11) 

0.5630(11) 2.06(23) 0.1797(12) 0.2247(12) 

0.2406(13) 3.22(24) 0.7007(13) 0.350413) 

0.1899(14) 2.58(25) 0.1113(15) 0.2785(14) 

0.5438(14) 2.06(26) 0.1729(16) 0.2818(15) 

0.8791(14) 3.22(27) 0.9657(17) 0.3184(16) 

0.1022(16) 4.42(27) 0.1667(18) 0.4212(16) 

0.1489(16) 5.15(27) 0.2180(18) 0.8701(16) 

0.2453(16) 5.24(27) 0.2241(18) 0.8897(16) 

0.2466(16) 5.25(27) 0.2248(18) 0.8916(16) 
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